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Abstract. We analyzed a large number of RXTE/PCA 
observations of Cyg X-1 in the hard and soft spectral 
states with total exposure time of « 190 and « 10 ksec 
respectively and time resolution better than « 250 fis in 
order to investigate its variability down to a few millisec- 
onds time scales. The model of modifications of the power 
density spectra due to the dead time effect was tested us- 
ing RXTE observations of an extremely bright source - 
Sco X~l. The results of these tests demonstrated, that al- 
though some problems still remain, for sources like Cyg 
X-1 (which is an order of magnitude less bright than Sco 
X-1) present knowledge of the instrument is sufficient to 
reliably analyze the power spectra in the kHz frequency 
range. 

In both spectral states of Cyg X-1 we detected sta- 
tistically significant variability up to ~150-300 Hz with 
a fractional rms amplitude above 100 Hz at the level of 
~ 2 — 3%. The power spectrum in the hard state shows a 
steepening at frequency of 40 — 80 Hz with the power 
law slope changing from ^ 1.7 to ~ 2.3 — 2.4. In the soft 
state the slope of power spectrum changes from ~ 1 to 
~ 2 at the frequency of ~ 15 — 20 Hz without any evi- 
dence of further steepening up to ^ 100 — 150 Hz. These 
break frequencies represent the highest characteristic fre- 
quencies detected in the power density spectrum of Cyg 
X-1 so far. 

Key words: Accretion, accretion disks - Instabil- 
ities - Stars:binaries:general ~ Starsxlassification - 
Starsineutron X-rays: general - X-rays: stars 



1. Introduction 

The shortest timescales in the variability of the com- 
pact objects always was a very interesting topic of X- 
ray astronomy. The efforts to measure the highest fre- 



quency in the variability of Cyg X-1 can be traced back to 



rocket experiments in 70-ties (e.g. Rothschild et al. 1974 
piles 1981 ). It was found that Cyg X-1 demonstrates 



statistically significant variability of the X-ray fiux at 
least up to frequencies of ~10 Hz. Since then Cyg 
X-1 was extensively observed by various satellite mis- 
sions searching for millisecond and submillisecond time 



scales variability (HEAO-1, see 


e.g. Meckins et al. 1984 


Wen et al. 1996]; EXOSAT, e.g. 


Belloni & Hasinger 1990 


GINGA, e.g. 


Miyamoto & Kitamoto 1989). At the present 



time th e Rossi X-ray T iming Explorer observatory 
(RXTE, iBradt et al. 1993|) is the most powerful experi- 



ment for the investigation of the short time scale X-ray 
fiux variations. It combines high telemetry rate (up to 512 
kbps), large effective area of the detectors (~6400 cm^ in 
the soft X-ray band), high time resolution (down to 1/isec) 
and comparatively low deadtime distortions. The compar- 
ison of the characteristics of several X-ray experiments can 



be found in e.g. Giles et al. 1998 



The shape of the power density spectrum (PDS) of 
Cyg X-1 in its low/hard spectral state is well established 
(in the range from - lO^^ to -50-100 Hz) since EXOSAT 
observations - it can be roughly represented by a constant 
from ~ 10~^ Hz up to some break frequency {P ^ /°), 
then it steepens to the power law with index ^1 {P 
f~^) and then it breaks again to the power law with the 
index -1.6-1.8, P ~ f-^ e-i s (gge e.g. RXTE results in 
Nowak et al. 1999|) . Short observations (—10 ksec) of Cyg 
X-1 by RXTE/PCA with high time resolution (- 4^s) 
showed the absence of significant v ariability of Cyg X -1 at 
time scales smaller than —3 msec ( Giles et al. 1998 ). 



In this paper we report the systematic study of very 
high frequency variability of Cyg X-1 in X-rays using large 
number of RXTE observations (total exposure time close 
to 190 ksec for the hard state and —10 ksec for the soft 
state) and carefully treating the deadtime effects. 
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2. Observations, data reduction and deadtime 
corrections 

For our analysis we used the data from the Proportional 
Counter Array (PCA) aboard RXTE. To study the source 
variability in the hard spectral state, we chose a compact 
group of observations of Cyg X-1 performed on 15-17 Dec. 
1996 (proposal P 10236) with total exposure time of ~80 
ksec and 56 available observations of Cyg X-1 from the 
proposal P30157 performed between Dec. 1997 and Dec. 
1998 with total exposure time of ssllO ksec. The soft spec- 
tral state data were accumulated on 4-18 June 1996 data 



(proposal P10512, see also Cui et al. 1997) with total ex- 
posure time of ~10 ksec. All observations were performed 
with 5 PCUs switched on. 

The power density spectra were calculated from the 
light curves with the bin duration equal to the intrin- 
sic time resolution of the data (2~^^sR::! 122/is in propos- 
als P10236 and P10512 and 2-'^'^s^ 244/is in proposal 
P30157). Some deviations from the noise level predicted 
by the PCA dead time models are observed in the power 
spectra constructed in the PCA energy sub-bands with 
the amplitude of ^ 10^® (rms/mean)^/Hz. We therefore 
followed recommendations of the RXTE GOF and ana- 
lyzed the light curves in the total PCA energy band, where 
these deviations seem to be considerably less significant 
(William Zhang, private communication). 

The power spectra calculated for individual time se- 
ries were normalized to units of squared fractional rms 



(Miyamoto et al. 1991) 



and averaged. In the above formula aj is the Fourier ampli- 
tude at the frequency fj ; N^- is total observed number of 
counts in the time series; -Rtot ~ averaged observed count 
rate for the time series. Note that the ideal Poissonian 
noise component is already subtracted, therefore Pj = 
for a Poissonian distribution of counts in the time series. 

In a real detector various effects could lead to modi- 
fication of the Poissonian noise level by an amount AP. 
Two major effects contribute to AP in the case of PCA 
detector: 



AP = AP 



dt 



AP.le 



where APjt is a modification of the noise level due to 
counts dead time, APvie - additional noise component re- 
sulting from vetoing of the PCA detectors by the Very 
Large Events (VLE). 

With an accuracy sufficient for our purpose the PCA 
dead time caused by the incident photons can be described 



as a non-paralizable process |^ ( [lahoda et al. 1996 ), for 
which modification of the noise level was derived by 
Vikhlinin et al. 1994 (see also ^hang et al. 1995 for the 
expression accounting for finite length of the time scries). 
We shall use below the Eq. (A4) from |Vikhhnin et al. 1994 
transformed to the units of squared fractional rms: 



APdt(rd, tb, P) = -4 • - . sin(2^/tb/2)2 x 
K 



k=k„ 



J2 X[2n{f + k/t^)]/{n{ft^ + k)f 



(1) 



k=-k„ 



where 



P'2[l - cos(27r/Td)] + P'27r/sin(27r/Td) 



R' 



P'2[l-cos(27r/rd)]2 
R 



[P'sin(2^/Td) + 2^/]2 



1-PTri 



Here R is the observed count rate, is the photon dead 
time, th is the bin duration, / is the frequency. For the 
practical purposes km of the order of 10 is sufficient. This 
equation (for km = oo) is valid for an infinitely long time 
sequence. One can use more complicated expression (see 



Eq. (44) from Zhang et al. 1995) if the analyzed time se- 
quence contains only small number of bins. 

To the first approximation the Very Large Events 
lead to appearance of an additional (positive) component 
(^hang et al. 1996]) 



APvlc(/, Tvlc, Pvlc) = 2PvleT^lc 



sin(7rTvie/) 



TTTvlc/ 



(2) 



where Tvie - the duration of the VLE window, Pvio - the 
VLE count rate in one PCU. Note that the above expres- 
sion does not include the binning and sampling effects and 
is valid in the limit PvieTVio << 1. 

Finally, the total noise component model in the first 
approximation would be: 



AP 



1 



[APdt (rd,ib,Pp cu/^vlc) /^vlc 



pcu 



1 



Mvlo 



1 Pvlc'^vlc 



where Ppcu - the total observed count rate in one 
pcu - number of PCUs. The APdt {Td,th,R) and 



PCU, iV, 

APvic (Tvie.Pvie) are given by the Eq. (1) and (2) respec- 
tively. The factor 1/A^pcu accounts for the fact that the 
count streams from A^pcu independent units were merged 

^ However, there is small paralizable co mponent in the total 
instrument deadtime ( Jahoda et al. 1996). The value of paral- 



izable deadtime in principle can slightly depend on the hard- 
ness of the source 
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Fig. 1. The power spectrum of Sco X-1 from the observa- 
tion performed with medium VLE window. The sohd 
hne shows the best fit model consisting of the noise level 
model AP and two kHz QPOs (see text). The dashed line 
shows the noise level component AP^t due to photon dead 
time. The lower panel shows the residuals data-model. 



together. The factor fi^ic accounts for reduction of the ob- 
served count rate due to Very Large Events. This formula 
is valid under the following conditions: i?vic << ^pcu, 

-Rvlc^vlo << 1 and Tvlc >> Tdt i?pcu (Tdt + l/i?pcu)- It 



is important to stress out that in the adopted dead time 
model i?pcu is the the total observed count rate of the 
events causing the dead time (as opposite to the ob- 
served count rate in the analyzed time series) . In the case 
of PCA detector it should include total Good Xenon rate, 
Propane and Remaining counts. We also note that for suf- 
ficiently small count rates, R ^3-5 kcnts/s/PCU depen- 
dence of APdt {Td,th, R) and, therefore, of AP upon R 
vanishes. The model has been tested on a series of Monte- 
Carlo simulations and has proven to be sufficiently accu- 
rate in the parameters range of interest. 

Note that we did not include in our model 
additional background term 

bright 



an 



Jernigan ct al. 200C ). For 



(Jahoda 1998 



sources 



(Jil 



or 
kc- 



nts/s/PCA) its contribution can be neglected. 

To verify our model for the noise component we com- 
pared it with the data of PCA observations of an ex- 
tremely bright source Sco X-1 (~100 000 cnts/s/PCA) 
for which all flavors of the dead time distortions are much 
more prominent than for Cygnus X-1. We used 3 differ- 
ent observations of Sco X-1 - 10059-01-01-00 (Feb. 14, 
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Fig. 2. The same as Figj^ but for observation #2, per- 
formed with short VLE window. 
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Fig. 3. The same as FigJ^ but for observation #3, per- 
formed with short VLE window and lower time resolution, 
w 122 IIS. 



1996; hereafter #1), 10059-01-03-00 (Feb. 19, 1996; #2) 
and 10056-01-02-02 (May 26, 1996; #3), having different 
value of the VLE window TvIc and different time resolu- 
tion. Ideally, once the values of the dead time Td and VLE 
window Tvlc are calibrated, the model should reproduce 
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Table 1. The expected and best fit model parameters of the deadtime correction for 3 observations of Sco X-1. 



Obs.ID. 






id, fJ-s 




i?pcu, kcnts/s 




Tvlc, MS 


-Rvlc 


, cnts/s 


x7dof 






exp. 


fit 


Tot^ 


GX= 


fit 


exp.'* 


fit 


exp.*" 


fit 




#1 
#2 

#3 


15.25879 
15.25879 
122.0703 


8-10" 
8-10" 
8-10" 


8.5 ±0.1 
8.7 ±0.1 
8.7(fixed) 


29.2 
19.6 
23.6 


20.8 
15.1 
18.3 


17.6 ±0.3 
10.0 ±0.5 
14 (fixed) 


150 
60 


152 ± 1 
76 ±2 
76.0(fixed) 


265 
239 
157 


274 ±2 
244 ± 13 
162 ±6 


152/151 
180/157 
43/46 





e.g. 



Jahoda et al. 199e, Jahoda et al. 1997, Jahoda 200C 



' - measured total count rate (Good Xenon, Propane Counts and Remaining Counts) per PCU 
- measured Good Xenon count rat e per PCU 



preflight values for all PCA, see Giles 1995 



measured count rate (HouseKeeping data) 



the noise level with count rates i?pcu and i?vie set to the 
measured values. However, as we discuss below, this is 
not strictly the case. We therefore followed the approach 



5x10-'' 



adopted by Jernigan et al. 2000 and fitted the power spec- 



tra at high frequencies, / > 300 Hz, with the model for 
the noise component plus the model for the source compo- 
nent leaving the noise model parameters free. The source 
contribution in this particular case was modeled as a su- 
perposition of two Lorentzians representing kHz QPOs. 
This approach has a disadvantage that it requires a pri- 
ory assumptions about the shape of the power spectrum 
of the source. In particular, if the source has very weak 
and very flat power spectrum component it might be not 
detected by our procedure. 

In Fig. H H and ^ we present the power spectra of Sco 
X-1 in observations #1 #2 and #3 along with the best fit 
noise models and residuals. The best fit parameters and 
their expected values are given in Table |l|. For observa- 
tion #3 we fixed parameters td, Rpcu and Tyie because the 
limited frequency range of the power spectrum (<4096 
Hz) does not allow us to constrain their values from the 
fit. From Fig. |l|, | and [s] one can see that with appro- 
priate tuning of the parameters the adopted noise model 
is capable to reproduce the observed shape of the noise 
component. 

In general the best fit parameters of the noise level 
model (Table |l|) look reasonable and are close to the ex- 
pected values. The best fit value of TvIc, in observation #2 
is by « 30% higher than the preflight value. It is possi- 
ble, however, that the duration of the VLE window differs 
from the preflight values (Alan Smale, private commu- 
nication). The most apparent difference is a large, by a 
factor of ~ 2, discrepancy between the best fit and ob- 
served total count rate. But we should note that the ob- 
served GoodXenon count rate is much closer to the best fit 
i?pcu value (see Table |l|). The reason of this discrepancy is 
not clear, it might be due to unaccounted processes in the 
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Fig. 4. The power spectrum of Cyg X-1 (observations 15- 
17 Dec, 1996; proposal P10236). Total exposure time ^ 
SOksec. Dashed line shows the instrumental noise level. 



PCA detectors at high count rates. We note however, that 
dependence of the noise component (in units of squared 
fractional rms) on the total count rate is a specific feature 
of a non-paralyzable dead time process at sufficiently high 
count rates. At lower count rate (i?pcu ^3000-5000 cnts/s) 
and in particular in the case of Cyg X-1 (i?pcu '--^1000-1500 
cnts/s) this dependence vanishes. 

3. Results 

We adopted the noise level model described in the previous 
section for the subsequent analysis of the Cyg X-1 data. 
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Fig. 5. The same as Fig. ||but for observations from pro- Fig. 6. The same as Fig. ^ but for observations from pro- 
posal P30157 (Dec. 1997-Dec. 1998, hard spectral state), posal P10512 (4-18 June, 1996, soft spectral state). Expo- 
Total exposure '^110 ksec. sure ~10 ksec. 



Since the Cyg X-1 observations were performed with time 
resolution coarser than w 122 fisec, the photon dead time 
and duration of the VLE window are not constrained by 
the data. We fixed them at the best fit values determined 
from the Sco X-1 data: = 8.7 fis and Tvie = 76/zs. As 
mentioned above, at the count rates typical for Cyg X-1, 



^1000-1500 cnts/s, the noise level does not depend 



on the count rate. We therefore fixed the i?pcu parameter 
at the averaged observed value. The only free parameter 
of the noise level model was i?vie- Similarly to the proce- 
dure applied to the Sco X-1 data, we included the source 
component into the model. The source component was 
represented by a power law P oc /~". We used the high 
frequency part of the power spectra, from 100 Hz up to 
the Nyquist frequency, for the fits. 

The high frequency part of the observed power spectra 
of Cyg X-1 before subtraction of the noise level AP{f), 
the best fit model and the residuals of the data from the 
model are shown in Fig. |4[ || and ^. The overall power 
spectra in a broad frequency range from 1 mHz to « 2 
IcHz after subtraction of the best fit noise level model 
are shown in Fig.^ One can see that the model describes 
the data points reasonably well: ~38/55 dof for obser- 



vations PI 0236 (Fig. 



■^56/49 dof for observations 



P30157 (Fig. |) and p ^46/49 for observations P1G512 
(soft state, Fig.^). The best fit values of the slope of the 
source power spectrum are: a — 2.4 ± 0.1 for P10236, 



a = 2.27 ± 0.08 for P30157 and a = 2.1 ± 0.2 for P10512. 
The intrinsic variability of the source has been statisti- 
cally significantly detected up to frequencies of ~150-300 
Hz (Fig. with the fractional rms amplitude in fre- 

quency band ^100-400 Hz of « 3% (hard spectral state) 
and «2% (soft state). 

In order to increase statistics we averaged all the hard 
state data increasing the total exposure time up to ~ 190 
ksec. The power law approximation of the averaged power 
spectrum above 100 Hz gives the power law index of 
a = 2.32 ± 0.07. The power law approximation of the 
same data in the 20-60 Hz frequency range results in 
a — 1.66 ± 0.01, although the PDS clearly has a more 
complicated shape than a power law (Fig.||). We therefore 
conclude that we detected statistically significant steepen- 
ing of the power spectrum of Cyg X-1 in the hard state at 
the frequency of ~40-80 Hz. Note, that some indication 
of such steepening can be found in Nowak et al. 1999| , al- 



though short duration of the used observations 20 ksec) 
was not sufficient to study the high frequencies in detail. 
In order to illustrate the high frequency behaviour of the 
power spectra, we plot in Fig.^ the power spectra multi- 
plied by p. As it can be seen from Fig.|| the particular 
shape of the rollover above 100 Hz is not significantly con- 
strained by the data. Assuming that power spectrum con- 
tinues with the same slope to the higher frequencies an 
observatory of the EXTRA/LASTE class (effective area 
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Fig. 7. The broad band (lO^'' — lO'' Hz) power spectrum 
of Cyg X-1 in the hard and soft spectral states averaged 
over different data sets used for the analysis. The power 
spectra were multiplied by the frequency. 



10 m^, Barret 2000 ) would be able to detect statisti- 



cally significant signal up to several kHz in ~tens of ksec 
exposure time and help us to say something more about 
the exact shape of the rollover of the PDS. 

The 2(7 upper limit on the amplitude of a possible QPO 
component (Lorentz profile with quality Q) in the 500- 
2000 Hz frequency range in the hard state power spectrum 
is ~2% for Q = 1 and ~0.9% for Q = 20. These upper Hm- 
its were obtained using the observations from the propsal 
P10236 performed on 15-17 Dec. 1996 in which the shape 
of the power density at the lower frequency did not vary 
significantly. The upper limit on the continuum noise com- 
ponent at the high frequency is somewhat more difficult to 
estimate. A statistical (2cr) error on the fractional rms in 
the 500-2000 Hz for the sum of all hard state observations 
is w 1%. Assuming that the noise level model is exact the 
upper limit on the fractional rms would coincide with the 
above number. However, as we mentioned in the previous 
section a weak and very flat source power spectrum com- 
ponent might be not detected given the procedure used to 
determine the instrumental noise level. 



4. Summary 

We analized ^190 ksec of observations of Cyg X-1 in the 
low/hard state and reanalyzed ~10 ksec of high/soft state 
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Fig. 8. The power density spectrum of Cyg X-1 at high 
frequencies. Note, that the power spectrum has been mul- 
tipHed by 

observations. We detected statistically significant variabil- 
ity up to frequencies of '^300-400 Hz with the fractional 
rms amplitude in the 100-400 Hz frequency range «3% 
in the hard state and r;2% in the soft state. The 2a up- 
per limits on the fractional rms amplitude of a QPO fea- 
ture with Lorentzian profile in the 500-2000 Hz frequency 
range are r^2% for Q - 1 and -0.9% for Q - 20 in the 
hard state and 2.5% (Q - 1) and 1.5% (Q - 20) in the 
soft state. Assuming that the model for the instrumental 
noise level is exact, the 2cr upper limits on the continuum 
variability in the 400-2000 Hz frequency range are —1% 
and — 1.7% for the hard and soft state respectively. 

We detected statistically significant steepening of the 
power spectrum in the hard state at frequencies '^40 — 80 
Hz. The power law approximation to the spectrum at 
frequencies higher than ^-^80-100 Hz gives the slope of 
a ~2.3-2.4. In the soft state the slope of power spectrum 
changes from ~ 1 to ~ 2 at the frequency of 15 — 20 
Hz without any evidence of further steepening up to the 
frequency of ^ 100 — 150 Hz. 

5. Astrophysical implications of high frequency 
variability 

Thus the sensitivity and time resolution of the RXTE 
allow one to probe interesting range of the variability time 
scales for accreting black holes. Several characteristic 
frequencies fall approximately in this range. First of all 
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the frequency of rotation at the last marginally stable 
orbit for the non-rotating (Schwarzschild) black hole is 



/ ~ 200 



IOMq 
M 



Hz. For the maximally rotating Kerr black 



hole this frequency ranges from '-^1200 Hz to ^ 150 Hz for 
the co-rotating and counter-rotating orbits respectively. 
If Keplerian rotation itself makes significant contribution 
to the variability of the X-ray flux (e.g. due to the 
Doppler shift and boosting of the emission from the "hot 
spots" in the disk) this difference should be manifested 
in this frequency range (Sunyaev 1972). Secondly typical 
sound crossing time of the disk height (in the case of 
optically thick geometrically thin standard accretion 
disk) corresponds to the frequency of hundreds of Hz. 
If strong turbulence is present in the accretion flow as 
anti cipated by the most popul ar theoretical models (see 
e.g. Balbusfc Papaloizou 1999 for review) then it may 
affect the variability at higher frequencies. The depen- 
dence of the amplitude of variations on the frequency 
is then a valuable tool to study the turbulence in the 
accretion flow. Thirdly light crossing time of the central 
region t\c = 3Rg/c ^ corresponds to the frequencies of 
the order of kHz. The light crossing time for the region of 
the main energy release (i? ~ 10-20 Rg) corresponds to 
the frequencies of 300-500 Hz. The character of the the 
variability may change drastically around this frequency 
or (e.g. in the comptonization models involving multiple 
scatterings of the photons) at a several times lower 



frequencies (e.g. Payne 1980, Sunyaev & Titarchuk 1980 



Pozdnyakov, Sobol, Sunyaev 198»^ 



Nowak fc Vaughan 199^ 

The RXTE data (i) demonstrate that power specifi- 
cally associated with any of these characteristic frequen- 
cies is not large (see given above upper limits on the nar- 
row features on the PDS in the frequency range from 
500-2000 Hz) - i.e. there are no "resonances" at these 
frequencies, (ii) limits the total power associated with 
high frequency variability (from 400 Hz to 2000 Hz) to 
a percent level, (iii) show that in the hard state of the 
source prominent, previously unknown, break in the PDS 
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